m 



PROCEEDINGS 

Of SCIENCE 



Non-equilibrium photon production arising from tlie 
chiral mass shift 



o 

(N 

u 

< 

00 






(N 
> 

m 

o 
o 

en 



% 



Frank Michler 

University of Frankfurt 

E-mail: piichler@th . physik . uni-f rankfurt . de 



Hendrik van Hees 

Frankfurt Institute for Advanced Studies 



E-mail: hees@f ias . uni-f rankfurt . de 



Dennis D. Dietrich 

University of Frankfurt 

E-mail: [dietrich@th . physik . uni-f rankfurt . de 



Stefan Leupoid 

Uppsala Universitet 



E-mail: Stefan . leupoldSphysics .uu.se 



Carsten Greiner* 

University of Frankfurt 



E-mail: carsten . greiner@th . physik . uni-f rankfurt . de 



We investigate the photon emission arising from the shift of the quark/antiquark masses during the 
chiral phase transition in the early stage of ultrarelativistic heavy-ion collisions. As this mass shift 
leads to spontaneous creation of quark-antiquark pairs and thus contributes to the formation of the 
quark-gluon plasma, our investigations are relevant in the context of finite-lifetime effects on the 
photon emission from the latter Earlier investigations on this topic by Boyanovsky et al. were 
plagued by a UV divergent contribution from the vacuum polarization and the non-integrability 
of the remaining contributions in the ultraviolet domain. In contrast to these investigations, we 
do not consider the photon number-density at finite times but for free asymptotic states obtained 
by an adiabatic switching of the electromagnetic interaction according to the Gell-Mann and Low 
theorem. This approach eliminates possible unphysical contributions from the vacuum polariza- 
tion and renders the resulting photon spectra integrable in the ultraviolet domain. It is emphasized 
that the consideration of free asymptotic states is indeed crucial to obtain such physically reason- 
able results. 
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1. Introduction and Motivation 



Ultrarelativistic heavy-ion collision experiments allow for studying strongly interacting matter 
under extreme conditions. One main objective of these experiments is the creation and exploration 
of the so-called quark-gluon plasma (QGP), a state of matter of deconfined quarks and gluons. The 
two basic features of the strong interaction, namely confinement ||I|] and asymptotic freedom |^ ^, 
predict that this state is created at high densities and temperatures [|[ ||, ||, ^, ||], which occur during 
ultrarelativistic heavy-ion collisions. 

The lifetime of the QGP created during a heavy-ion collision is in the order of up to 10 fm/c. 
After that, it transforms into a gas of hadrons. Thus, experiments cannot access the QGP directly, 
and the investigations of its properties have to rely on experimental signatures. One important 
category of these signatures are direct photons as electromagnetic probes. As they only interact 
electromagnetically, their mean free path is much larger than the spatial extension of the QGP. 
Therefore, they leave the QGP almost undisturbed once they have been produced and thus provide 
direct insight into the early stage of the collision. 

One important aspect in this context is that the QGP, as it occurs in a heavy-ion collision, 
is not a static medium. It first thermalizes over a finite interval of time, then keeps expanding 
and cooling down before it hadronizes finally. This non-equilibrium dynamics has always been 
a major motivation for investigations on non-equilibrium quantum field theory [^ |I0[ ^, 12, 13 



14, 15, 16, 17, 18, 19, 2C]. Besides the role of possible memory effects during the time evolution 
[0 m, H, |4|, H, H, Ivj, H, ||], it is of particular interest how the finite lifetime QGP itself 



affects the resulting photon spectra. 



The first investigations on this topic by Boyanovsky et al. [30, 31, B2] have shown that the 



finite lifetime of the QGP gives rise to the contribution to the photon yield from first-order QED 
processes, i.e., processes of linear order in the electromagnetic coupling constant, a^, which are 



kinematically forbidden in thermal equilibrium. Furthermore, the results from [ )30| , pl| , [32| ] sug- 
gested that such non-equilibrium contributions dominate over leading-order thermal contributions 
in the ultraviolet (UV) domain. The latter contributions ai^e linear both in the electromagnetic 
coupling constant, a^, and the strong coupling constant, a^, and hence of overall second order. 



On the other hand, however, the investigations in [30, BJl, 32] had been accompanied by two 



artifacts. In the first place, the vacuum polarization was found to lead to a divergent contribution 
to the photon number-density for given photon energy. Furthermore, the remaining contributions 
to this quantity did not decay fast enough with increasing photon energy for being integrable in 
the UV domain. In particular; the total number density and the total energy density of the emitted 
photons (after subtracting the contribution from the vacuum polarization) were logarithmically and 
linearly divergent, respectively. 

Later on, the topic was also picked up by Fraga et al. [|3|, 34 ] where the ansatz used in 



31, 32] was considered as doubtful, as it came along with the mentioned problems. In particular. 



the concerns raised in [ )33| , [34| ] were the following: 

] the time at which the photons are observed has been kept finite. Either this 



31 



In[ 

corresponds to measuring photons that are not free asymptotic states or it corresponds to 

suddenly turning off the electromagnetic interaction at this point in time. Both cases are 

questionable. 
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The initial state in [30, 31, B3] has been specified for a system of themialized quarks, an- 



tiquaries and gluons not containing any photons. On the other hand, however, a system of 
quai^ks and gluons, which undergoes electromagnetic interactions, necessarily contains pho- 
tons. Hence, taking an initial state without any photons and without the Hamiltonian for 
electromagnetism corresponds to switching on the electromagnetic interaction at the initial 
time, which is questionable as well. It was shown in |35] that the ansatz used in [|^, 31, 32] 
is indeed equivalent to such a scenario. 



The divergent contribution from the vacuum polarization is unphysical and thus needs to be 



renormalized. Nevertheless, the renormalization procedure presented in [32] is not coherent 
since no derivation of the photon yield with rescaled field operators has been presented in 



[32]. 



The authors of [ p3[ [34| ] did, however, not provide an alternative approach for how to handle the 



mentioned problems in a consistent manner. Solely in [ p5[ ] it was indicated that the question of 
finite-Ufetime effects could be addressed within the 2PI (two-particle irreducible) approach even 
though the conservation of gauge invariance remains challenging. Later on Boyanovsky et al. 



insisted on their approach [36] and objected to the ai^guments by [33, 34] as follows: 



Non-equilibrium quantum field theory is an initial-value problem. This means that the den- 
sity matrix of the considered system is first specified at some initial time and then propagated 
to a later time by the time-evolution operator. For that reason, the Hamiltonian is not modified 
by introducing a time-dependent artificial coupling as it would be the case for a 'switching 
on' and a later 'switching off of the electromagnetic interaction. 

The quark-gluon plasma, as it occurs in a heavy-ion collision, has a lifetime of only a few 
fm/c. Therefore, taking the time to infinity is unphysical as this limit requires the inclusion 
of non-perturbative phase-transition effects on the photon production. 



The renormalization technique of [ pq ] provides a rescaling of the photon field operators such 
that the photon number operator actually counts asymptotic photon states with amplitude 
one. 



During this debate, however, the original problems encountered in [ |30[ , |31| , |32[ ] remained open. 
In order to resolve them in a satisfactory manner, we have followed two different approaches. 

In the first approach ^%, we have modeled the finite lifetime of the (thermalized) quark-gluon 
plasma during a heavy-ion collision by time-dependent occupation numbers in the photon self- 
energy. By means of this procedure, we have been able to renormalize the divergent contribution 
from the vacuum polarization in a consistent manner. However, it has not been possible to get the 
problem with the UV behavior fully under control. 

We have expected in the first place that this remaining shortcoming results from a violation of 
the Waid-Takahashi identities. For that reason, we have also considered a conceptionally different 



scenario [38] where the production of quark-antiquark pairs together with photons results from 
a change in the quark-antiquark mass. Such a change occurs during the chiral phase transition 
in the very early stage of a heavy-ion collision. It has been shown in [^, ^] that it leads to 
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the spontaneous and non-perturbative pair production of quarks and antiquarks, which effectively 
contributes to the formation of the QGP. We have investigated the photon production arising from 
this pair-creation process. As this photon production is effectively induced by the change of the 
quark/antiquark mass and thus by the chiral phase transition, it is from now on referred to as chiral 
photon production. 



In contrast to [37|, such a scenario has the crucial advantage that it allows for a first-principle 
description by introducing a Yukawa-like source term in the QED Lagrangian. This source term 
couples the quarks and antiquarks to a time-dependent, scalar background field. In this way, they 



obtain a time-dependent mass, which conserves the Ward-Takahashi identities. As in [|37|], we have 
restricted ourselves to first-order QED processes. They are kinematically possible since the quarks 
and antiquarks obtain additional energy by the coupling to the background field. To the contrary, 
the coupling to the background field is resummed to all orders as to properly take into account the 
non-perturbative nature of the pair-creation process. 



In this context, there is another crucial difference to the approaches in [gO[ |31|, ^ ^7|] : There 
the photon number-density has been considered at finite times. In the course of our investigations 
[^8|], however, we have shown that the photon number-density has to be extracted in the limit f — > oo 
for free asymptotic states, which are the only observable ones since they reach the detectors. Such 
states have been obtained by introducing an adiabatic switching of the electromagnetic interaction 
according to the Gell-Mann and Low theorem. The photon number-density is then considered 
in the asymptotic time limit and the switching parameter is taken to zero at the very end of our 
calculation. 

We shall demonstrate that this procedure eliminates a possible unphysical vacuum contribution 
and, furtheiTnore, leads to photon spectra being integrable in the ultraviolet domain if the time 
evolution of the quark/antiquark mass is modeled in a physically reasonable manner. It is also 
emphasized that the consideration of the photon number-density for free asymptotic states is indeed 
essential to obtain such physically reasonable results. 

2. Model description 

Before we turn to our first-principle description on chiral photon production, we first provide 



a short insertion on our model description [[37||, where we have made an ansatz on the two-time 
dependence of the photon self-energy. For this purpose, we have taken into account that the vac- 
uum contribution to the photon self-energy is always persistent, whereas the medium contribu- 
tion only occurs as long as the QGP is actually present. These two aspects had been disregarded 



in pOl , |31| , |32[ ]. This time dependence has been in implemented by introducing time-dependent 



quark/antiquark occupation numbers 

nF{E)^nF{E,t)=f{t)nF{E), (2.1) 

in the photon self-energy. Here nY^{E) denotes the Fermi-Dirac distribution at the given tempera- 
ture, T, and f{t) is the function modeling its time evolution. This function changes monotonously 
from zero for t -^ — oo (vacuum) to one for ? — )• oo (QGP fully persistent) over an assumed formation 



time interval, T. As in [3C, 31 , 32], we have considered the photon number-density at finite times. 
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By means of this procedure, we have been able to renormalize the divergent contribution from 
the vacuum polarization in a consistent manner. For our numerical investigations on the remain- 
ing medium contributions, we have compared the resulting photon spectra for different switching 
functions, /,(?), which are depicted in Figure |I} As one can see, /i(f) describes an instantaneous 
formation of the QGP at f = 0, whereas /2(0 describes a formation over a finite interval of time, t, 
for which we have assumed T = 1 .0 fm/c. 







f,(l) — 
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Figure 1: The time evolution of the QGP is modeled by different switching functions, fi{t). 

It has been our initial hope that in addition to the renormalization of the vacuum contribution 
our model description would also lead to UV integrable photon spectra if the the formation of the 
QGP is assumed to occur over a finite interval of time. Such an assumption is reasonable from the 
phenomenological point of view. A comparison of the photon spectra for the different switching 
functions, which is done in Figure ^ does, however, show that this is only partly the case. 




10 12 14 16 18 20 
tOk [GeV] 



Figure 2: Comparison of the photon spectra for different switching functions. For /i(f) describing an 
instantaneous formation, the photon spectrum decays as l/o^ and exhibits only a slightly steeper decay 

ex 1/(0^-^ for flit)- Or ~ \k\ denotes the photon energy for given photon three-momentum, k. 



For the case of an instantaneous formation the photon number-density scales as l/ft)^ (with 
ft)^ = 1^1 denoting the photon energy for given photon three-momentum, k) in the UV domain. 
This implies that the total number density and the total energy density of the radiated photons are 
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logarithmically and linearly divergent, respectively. This artifact is only partly resolved if we turn 
to a fonnation over a finite interval of time describing a physically more realistic scenario. In this 
case, the UV scaling behavior of the photon number-density is only suppressed to 1/ft)?'^ such that 
only the total number density is rendered UV finite whereas the total energy density remains UV 
divergent. For our numerical investigations, we have chosen /2(f) to be continuously differentiable 
once but one can show that this artifact persists for arbitrarily smooth switching functions. 

We have assumed in the first place that the still problematic UV behavior arises from a viola- 
tion of the Wai^d-Takahashi identities for the photon self-energy within our model description [^]. 
Therefore, we aimed to find a new approach in which these identities are conserved. For that rea- 
son, we have eventually studied chiral photon production as this scenario allows for an accordant 
first-principle description. 

3. Chiral photon production 

We model the change of the quark/antiquark mass during the chiral phase transition by intro- 
ducing a Yukawa-like source temi in the QED Lagrangian 

^{x)=^QEB{x)-g(P{t)Y{xmx) . (3.1) 

The source field, <p{t), is assumed to be classical and time-dependent only, which effectively assigns 
the quarks and antiquark a time-dependent mass 

m{t)=mc+g(p{t) . (3.2) 

Since such an ansatz keeps the Lagrangian gauge invariant it conserves the Ward-Takahashi iden- 
tities of QED, which we will demonstrate explicitly below for the photon self-energy. 

The change of the quark/antiquark mass from its constituent value, m^, to its bare value, mt 
leads to the spontaneous and non-perturbative pair creation of quarks and antiquarks. We investi- 
gate the photon emission resulting from this pair-creation process. Thereby, we assume that our 
system does not contain any quarks/antiquarks or photons initially. The initial density matrix is 
hence given by 

p{to) = \Oyci)^\Oy) . (3.3) 

Here \Oqq) and |0y) denote the vacuum states of the fermionic and the photonic sector, respectively. 
As the quark/antiquark mass changes with time, it is important to point out that the former is 
defined with respect to the initial constituent mass, m^. The initial time, to, is chosen from the 
domain where the quark/antiquark mass is still at this value, i.e., tQ < t^ with t^ denoting the time 
at which the change of the quark/antiquark mass begins. In the case of parameterizations, m{t), 
where the time derivative, rh{t), has a non-compact support, it is sufficient to ensure that one has 
|g(/)(0l<mcfor?<4 

Since the quark mass is time-dependent only and our system is initially given by the vacuum 
state (^J|), it is spatially homogeneous. For such a system, the photon number-density is formally 
given by 
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The sum runs over all physical (transverse) polarizations and the average is taken with respect to the 



initial density matrix ([3.3|). Moreover, d{k,X,t) together with its Hermitian conjugate corresponds 
to the expansion coefficients in the plane-wave decomposition of the photon-field operator, A^ {x,t), 
i.e., 




e^{k,X)a{k,X,ty-' + s*{k,X)d'^{k,X,t)e 



^t/ 



-ikx 



(3.5) 



Expression ( |3.4D for the photon number-density has also been used in pOj , pT] , p2| ] with the initial 
density matrix, p(fo). being instead specified for a thermalized QGP not containing any photons. 



However, before one can start with any further evaluation and numerical analysis of (3.4), one 
has to clarify under which circumstances d(k,X,t) and its Hennitian conjugate actually allow for 
an interpretation as a single-photon operator. Such an interpretation is not given in general since 
A^{x,t) describes an interacting electromagnetic field. It is, however, possible in the asymptotic 
time limit ? — )• oo for free asymptotic fields. Such fields ai^e obtained by introducing an adiabatic 
switching of the electromagnetic interaction according to the Gell-Mann and Low Theorem, i.e., 

HEM{t) = eJd\p{x)Y^Y{x)A^'{x)^e-'\'\HEM{t), with £>0. (3.6) 

In order to obtain physically well-defined results for the photon number-density, we have to specify 
our initial state for ?o -^ —°° and consider (3^) in the limit ? — ;• oo for free asymptotic states. 
This procedure is illustrated in Figure ^ Since the introduction of an adiabatic switching of the 
electromagnetic interaction per se is an artificial procedure we have to take e — >• at the very end 
of our calculation. 




Figure 3: Schematic representation of our asymptotic description. The initial state is specified for to — ^ —°° 
(a). From the adiabatic switching of the electromagnetic interaction, one then has interacting fields around 
f = (b) which evolve into free fields in the asymptotic time limit r — > oo (c), where the photon number- 
density is considered. Only at the very end of our calculation, the switching parameter, e, is taken to zero. 



As the electromagnetic coupling is small, we evaluate ( ^^ to first order in a^, but keep all 
orders in g. The latter is achieved by constructing an interaction-picture representation that incor- 



porates the source term [p8|]. Together with (3.6), this procedure leads to 

d?i /' df2in^(^,fi,?2)e''"*'^""'^^ 



d^nUt) 
Itth , \ 



1 



\Z.7L] J — oo J — oo 



(3.7) 
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with the underUning denoting that 



du= dwe^'^l"! . 

D J — OO 

Moreover, illj(k,ti,t2) describes the transverse part of the photon self-energy, i.e., 

in<{kA,t2) = f\k)in<^{k,h,t2) ■ (3.8) 

Y^^{k) is the photon polarization tensor reading 

f-T]^^-^ for jU,VG {1,2,3} 
Y^'ik) = £ e^'-*ik,X)s'{k,X) = \ ''i ' \ (3.9) 

A=± I Otherwise 

with Tj'^^ = diag{l, — 1, — 1, — 1}. The photon self-energy, illl2y(k,ti,t2), is in turn given by the 
one-loop approximation 

in<^{kAJ2) = e^ J T^Tr{^y^S^{p + k,tut2)YvSUp,t2,h)^ , (3.10) 

with the propagators fulfilling the equations of motion 

[ifd,,+Ypi-mih)]sfipA,t2) =0, (3.11a) 

[ifd,^-Ypi + m{t2)] sfipA,t2) = . (3.11b) 

It follows from these equations that the photon self-energy ( p.lO| ) fulfills the Ward-Takahashi iden- 
tities 

4, iU^^ (k, HA)- i^^'in<, (^, f 1 , f2) = , (3. 12a) 

d,^U<Q{kA , t2) + ikJin<j{k, tut2) = 0. (3. 12b) 

The one-loop approximation for the photon self-energy includes the processes of first order in 
«£, which is shown in Figure ^. In particular, these processes are (one-body) quark/antiquark 
Bremsstrahlung, quark-antiquark pair annihilation into a single photon and the spontaneous cre- 
ation of a quark-antiquark pair together with a photon out of the vacuum. They are kinematically 
possible since the quarks and antiquarks obtain additional energy by the coupling to the background 
field, (/)(?)■ 



As the above diagrammatic arguments indicate, the photon number-density (3.7 ) can be written 



as the absolute square of a first-order QED transition amplitude and, as a consequence, is positive 



semidefinite. To demonstrate this explicitly, we first undo the contraction (3.8) in (3.7), which leads 
to 



2«^rof = 7^ ^ d^in<,(^,^i,?2)e-^("-^). (3.13) 



As the next step, we expand the fermion propagators in terms of positive- and negative-energy 
wavefunctions, V'iJ.itilO for given momentum, p, and spin, s, i.e., 

S^{pA,t2) = -i£'/^p,.,t(^i)m.v,t(f2) , (3.14a) 

.V 

S^{pA,t2) = i£'/^v,;(fi)'/>p,a(^2) • (3.14b) 
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Figure 4: The photon self-energy, ill.^ {k, fi , r2 ), is given by the one-loop approximation including the first- 
order QED processes. These processes are quark Bremsstrahlung (a), antiquark Bremsstrahlung (b), quark- 
antiquark pair annihilation into a single photon (c), and the spontaneous creation of a quark-antiquark pair 
together with a photon out of the vacuum (d). These processes are possible because of the coupling of the 
quarks and antiquarks to the background field, (j){t). This coupling is depicted by the open circles as to 
distinguish it from the electromagnetic coupling depicted by the full circles. 



These wavefunctions fulfill the equation of motion 

[ifdt + Ypi - m{t)] Yp,s;nit) , 
with the initial conditions 

V^P,.,t(0 ^ V^,.,t(0 for t<t'o, 
Vp^it) ^ WpM^t) for ?<?o. 



(3.15) 

(3.16a) 
(3.16b) 



Here V^^fi (0 denote the positive- and negative-energy spinors with respect to the initial, con- 



stituent mass, jjic. Upon insertion of ( [3.8[ ) and relations ( |3.14| ) into ( |3.13| ), we can finally rewrite 
the photon number-density as 

2 



2C0r- 



* d^xd^k 






d^p 



due^{k,X)- Yp^,.d'')rti¥p+ksi("y"^'" 



(3.17) 



This absolute-square representation ensures that (3.7) cannot acquire unphysical negative values 



The physical photon number-density is extracted from (3.7) by taking the successive limits f — > oo 
and e — )• after the time integrations have been peiformed, i.e.. 



2(Ot 



d^n 



r 



''d^d^k 



lim 



1 



(2;r)3 7-00 — 7-00 — 



d?i / dt2iUj{k,ti,t2)e 



iWt{ti-t2) 



(3.18) 



in^(^,fi,?2) reduces to the vacuum polarization if both time arguments are taken from the domain 
where the quark/antiquark mass is still at its initial constituent value, rric, 



2e' 



d^p 

(2^ 



1 + 



px{px + k)-\-m^ 

P p+k 



^ leH^^+*+^^J("-'^) 



(3.19) 
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From the mass-shift effects, in^(^,?i,?2) will acquire an additional non-stationary contribution, 
i.e., 

in<(k,ti,t2) = in<o(kA-t2) + iAn<{k,ti,t2) , (3.20) 

depending on both time arguments explicitly. When performing the successive limits f — )• oo and 



£ — > the contribution from the vacuum polarization (3.19) is eliminated. This can be seen by 



inserting (3.19) into (3.18), which leads to 



(Of. 



d^ny 



kA'i 



d^xd^k 



VAC 



:lim--— -r 

<lim--— -r 

e^o {Inf 



(2;r)3 






1 + 



1 + 



px{px + COj^) + ml 

P p+k 



2e 



px[px + (Of) + m\ 



P p+k 



e'^+[Ei ^ + El + (Of 

p+k P k 



£^ -. + Ei + (Of 

p+k P k 



=0, 



(3.21) 



where we have taken into account that Ei -.+£■£ + ft)^ is positive definite in the second step. There- 



fore, only mass-shift contributions to ( 3.18 ) characterized by iAITj (^, ?i , ?2 ) remain. In this context, 
we would like to point out that adhering to the exact sequence of limits, i.e., taking first f — > oo and 
then £ — )• 0, is indeed crucial to eliminate the contribution from the vacuum polarization. When 



first taking £ ^ at some finite time and f — )• oo afterwards, the contribution from (3.18) instead 
turns into 

px{px + COf) + ml I 



d^ny(t) 

(Or ^ '\ 



VAC 



{iTiy J (271 



dV 



1 + ' 



1 



E'iE'^ -. 

P p+k 



Ei ^ + Ei + (Of 

p+k P k 



2 ' 



(3.22) 



with the integral over d^p being linearly divergent for given photon energy, (O^. In the course of 
our numerical investigations, we furthermore demonstrate that the conect sequence of limits is also 



essential to obtain physically reasonable results from the mass-shift effects. Together with (3.10) 



and ( p.lll ), expression ( |3.18| ) describes photon emission arising from the chiral mass shift at first 
order in a^ but to all orders in g. 



4. Quark-antiquark pair production 



Before we turn to our numerical investigations on chiral photon production, we first provide 
an insertion on quark-antiquark pair production arising from the chiral mass shift. Thereby, we 
consider the pair production arising due to the chiral mass shift only. The starting point of our 
investigations hence is the fermionic part of the Hamiltonian 



H{t)= d^xp{x) -iY-V + m{t) \if{x) , 



(4.1) 



with the mass function, m{t), given by ( |3.2| ). As the next step, we expand the fermion-field operator, 
\j/{x), in terms of the positive- and negative-energy wavef unctions, YpAli'')' ^•^•' 



Vi^) = L 






¥p,s,t{t)h,^ + ¥pM(t)^- 



p,s 



^ip-x 



(4.2) 
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with bps and d^ps both annihilating the initial fermionic vacuum state, \0q^). Upon insertion of 
(^) into dOl ) we obtain 



^(0 = l/ (2^ {^(^) hshs - d-p,dl.,,] +A{t)bl/_.^^^ + A*{t)d\p/bp,} . (4.3) 
In order to keep the notation short, we have introduced 

a{t) = V>p,.. t(0 [Ypi + m{t)] V//j,,,t(0 = -¥pMi^) [YPi + m{t)] Yp^i^) > (4-4a) 

A(0 = '^^..tlO [YPi + m{t)] YpMit) ■ (4-4b) 

The particle number-density of quarks and antiquarks, which coiTcsponds to the number of pro- 



duced quark-antiquark pairs per unit volume, is extracted from ( [4 .31 ) by diagonalizing this expres- 
sion via a Bogolyubov transfonnation 



n 



bp,s{t) = ^pAt)bp,s + Tlp,,{t)dl^^^ 



dK^it) = ^^t)dl.^^-r^lMb 



p.s 



(4.5a) 
(4.5b) 



As discussed in greater detail in [i38|], this procedure corresponds to a re-expansion of the fermion- 



field operator ( [4.2[ ) in terms of the instantaneous eigenstates of the Hamilton density operator 

hv>{t) = -iy-y + m{t) . 
As a consequence, the Hamiltonian ( |4. 1[ ) reads in terms of the transformed creation and annihilation 



operators ( \.5 ) 



dV 



^(0 =1/ (43^.^(0 [bW^pAt)-d-pAt)dl.,{t) 



dV 



Ep{t) bl {t)bpAt) + dK {t)d^pAt) 



(4.6) 



In the second step, H{t) has been normal-ordered with respect to ( |4.5| ) as to avoid an infinite 
negative vacuum energy. Moreover, we have introduced the dispersion relation 



Ep{t) 



-nfi{t) . 



(4.7) 



The Bogolyubov particle number-density and thus the number of produced quark/antiquark pairs 
is then defined as pT] ] : 



■LM^~^-pAf)^^pAt)K) 



1 



{2n) 



Ll%^-(0| 



(4.8) 



For our numerical investigations on the number of produced quark-antiquai^k pairs, we consider 
different mass parameterizations, which are depicted in Figure ^ 
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Figure 5: The change of the quark/antiquark mass from its constituent value, m^, to its bare value, mt,, 
during the chiral phase transition is modeled by different mass parameterizations, m,(f ). 



Thereby, mi{t) describes an instantaneous change from the constituent mass, m^, to the bare 
mass, nih, at f = 0, whereas m2{t) and mT,{t) each describe a change over a finite interval of time. 
The latter two mass parameterizations differ with respect to their order of differentiability, mjit) 
is continuously differentiable once, whereas wi3(f) is continuously differentiable infinitely many 
times. For both parameterizations, we assume a transition time of T = 1 fm/c. Figure ^ compares 
the asymptotic particle spectra for the different mass parameterizations, m;(?). 




0.2 0.4 0.6 0.8 I 1.2 1.4 1.6 l.S 

p [GeV] 



Figure 6: The decay behavior of the asymptotic particle spectra shows a strong sensitivity on the 'smooth- 
ness' of considered mass parametrization, m,(f). 

For mi{t), the quark/antiquark occupation numbers scale as l/p^ for p ^ mc,mh, which im- 
plies that the total number density and the total energy density of the produced quarks and anti- 
quarks are linearly and quadratically divergent, respectively. However, this artifact is cured when 
turning from instantaneous mass change to a mass change over a finite interval of time. In partic- 
ular, the quai^k/antiquai^k occupation numbers are suppressed to l/p^ for niiit) such that both the 
total number density and the total energy density are rendered finite. Moreover, if one turns from 
m2{t) to m3(?), which is continuously differentiable infinitely many times and hence describes the 
most physical scenario, the quark/antiquai^k occupation numbers ai^e suppressed even further to an 
exponential decay. One encounters the same sensitivity on the 'smoothness' of the mass parame- 
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terization if the mass is not only changed from its constituent value, nic, to its bare value, m/„ but 
also back to its constituent value after a certain interval of time as to take into account the finite 



lifetime of the chirally restored phase [|38|]. 

In addition to the asymptotic particle spectra, we also consider the time evolution of the 
quark/antiquark occupation numbers, which is shown in Figure ^ 



p=0.5GeV / 


\ ' 


10 






p=I.OGeV / ^ 




10- 


p=1.5GeV / /\ 






p=2.0GeV / //\ 


r\ _ 




\Aj^ - 


10 


M^- - 


^ 1^-1 


n yvN^-- ■ 


'^ 10"' 


//// 




10' 


/// m(t)^m 2(1) 


T=1.0fm/c 




/ m =0.35 GeV 




10"' 


' m^sO.OI GeV 


1 


I0-' 




t [fm/c] 



t [fni/c] 



Figure 7: Time dependence of the quark/antiquark occupation numbers for m^it) (left panel) and m^(t) 
(right panel) for T = 1 .0 fm/c. For large values of p, they exhibit an 'overshoot' over their asymptotic value 
around f = 0, which is particularly distinctive for mT,(t). 



As one can see, the quark/antiquark occupation numbers exhibit a strong 'overshoot' over their 
asymptotic value in the region of strong mass gradients, which is especially distinctive for m-iit). 
In particular, the quark/antiquark occupation numbers scale as Xjp^ for large p when the mass 
is being changed from mc to m\j. Such a scaling behavior means that the total number density is 
finite, whereas the total energy density is divergent. This divergence, however, only shows up in 
the region of strong mass gradients and disappeai^s again as soon as the quaiVantiquark mass has 
reached its final bare value, m/,. 

At first sight, the temporary logarithmic divergence might be disturbing. In this context, it is, 
however, important to point out that only the asymptotic energy density, i.e., for ? ^ 00 constitutes 



a well-defined physical quantity since the interpretation of ( [4.8[ ) as a particle number-density is 



only justified for asymptotic times where m(f) = 0. The reason is that the dispersion relation (^.7) 
then actually characterizes free and thus detectable particles, wheres it describes quasiparticles for 
m{t) 7^ 0. In particular, the asymptotic energy density is finite as long as the mass parameterization, 
m{t), is chosen smooth enough, which represents a physically reasonable condition. Moreover, the 
temporary logarithmic divergence in the energy density does not manifest itself in form of a similar 
pathology neither in the total number density nor in the total energy density of the emitted photons. 



5. Photon production 

Our investigations on pair production have shown that the asymptotic quark/antiquark occu- 
pation numbers exhibit a very strong sensitivity on the 'smoothness' of the mass pai^ameterization, 
m{t). In particular, the asymptotic particle spectra are rendered UV integrable if the mass change 
is assumed to take place over a finite interval of time, T. We now determine whether the asymptotic 
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photon number-density exhibits a similai^ sensitivity. In this context, we recall that only asymp- 
totic quantities constitute observables whereas the corresponding expressions at finite times do, in 
general, not allow for a similar interpretation. 

For our investigations on chiral photon production, we again consider the different mass pa- 
rameterizations shown in Figure |[ For the special case of an instantaneous mass shift, the loop 



integral entering ( [3.18[ ) features a linear divergence. This divergence arises from the Bogolyubov 
particle numbers oc l/p^ for large p [38]. This behavior is an artifact from the instantaneous mass 



shift and is removed if the mass shift is assumed to take place over a finite interval of time. Hence, 
from the conceptional point of view this divergence does not require a specific kind of renormal- 
ization and is regulated by a cutoff at p = Ac. Figure 7 shows the resulting photon spectra for the 
different mass parameterizations. Thereby, we have chosen Ac = 20 GeV for m{t) = mi{t). 
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Figure 8: The decay behavior of the asymptotic photon spectra in the UV domain is highly sensitive to the 
'smoothness' of the mass parameterization. In particular, they are rendered integrable in this domain if the 
change of the quark/antiquark mass is assumed to take place over a finite interval of time, T. 



We see that for an instantaneous mass change, the photon number-density scales as I/oj? for 
large (Oi^, which means that for given given Ac, the total number density and the total energy den- 
sity of the emitted photons are logarithmically and linearly divergent, respectively. In contrast to 



our earlier model approach [37], this artifact is now fully removed if we turn from an instanta- 
neous mass shift to a mass shift over a finite interval of time. For rn2{t), which is continuously 
differentiable once, the photon number-density is suppressed from l/(o} to l/o^ in the UV do- 
main such that both the total number density and the total energy density of the emitted photons 
are UV finite. Moreover, if one turns from ni2{t) to rnT,{t), which continuously differentiable in- 
finitely many times, the photon number-density is suppressed even further to an exponential decay. 
The sensitivity of this quantity to the 'smoothness' of the mass pai^ameterization is the same if the 
quark/antiquark mass is changed back to its constituent value, Mc, after a certain period of time as 
to take into account the finite lifetime of the chirally restored phase [|^. 



As the artifact with the UV behavior encountered in [30, Bll, 32] and still partly in [37] is 



now resolved, it is convenient to compare our results to leading-order thermal contributions. This 
is done in Figure ^ Thereby, the thermal contributions have been obtained by integrating the 



leading-order thermal rate taken from [i2\ over the expected lifetime of the chirally restored phase 
of Tl = 4.0 fm/c at a temperature of T = 0.2 GeV. In this context, we note again that first-order 
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photon production does not occur in thermal equihbrium. There the leading-order contributions 
are of linear order both in the electromagnetic coupling constant, a^, and the the strong coupling 
constant a^, and hence of overall second order. 




Figure 9: Comparison of first-order mass-shift contributions to leading-order thermal contributions. For 
/ = 3 most likely characterizing a physical scenario, the contributions arising from the chiral mass shift are 
clearly subdominant for ft^ ^ 1 GeV. 



Nevertheless, we see that chiral photon production does not generally dominate over leading- 
order thermal photon production. For m{t) = mj^t), which describes the most physical scenario as 
it is continuously differentiable infinitely many times, we see that chiral photon production is even 
subdominant by several orders of magnitude in the UV domain. 

For completeness, we would still like to demonstrate that adhering to the correct sequence 

oo and then e 



of limits leading to ( |3.18| ), i.e., taking ^jr^f ? — )• oo and then £ — )• is crucial not only to eliminate 
a possible unphysical contribution from the vacuum polarization, but also to obtain physically 
reasonable results from mass-shift effects. For this purpose, we consider the time evolution of the 



pure mass-shift contribution to (3.18) for different values of e which is shown in Figure 10 





Figure 10: Time evolution of the mass-shift contribution for different values of e for m2{t) (left panel) and 
mj{t) (right panel). 

Taking first £ — )• at some finite time corresponds to the cui^ve labeled by £ = in each case. 
If we then take the subsequent limit ? — > oo, we obtain an asymptotic value for the mass-shift contri- 
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bution which is by several orders of magnitudes larger than the value against which the asymptotic 
values for finite e converge in the limit £ — )• and which accordingly corresponds to ( ^.18 ). The 



reason why the asymptotic value obtained for an interchanged sequence of limits must be unphys- 
ical is that it coincides with the transient value of the mass-shift contribution for sufficiently small 



values of e. On the other hand, expression (3.7) does not allow for an interpretation as a photon 



number-density at finite times since one either has no free asymptotic states or one has to intro- 
duce an artificial switching-off of the electromagnetic interaction at the point of time at which 



(3.7) is considered. By means of the latter procedure, a violation of the Wai^d-Takahashi identities 



furthermore comes into play again [38| 



The same conceptual problem arises when only using an adiabatic switching-on of the elec- 
tromagnetic interaction for t — )■ — oo but no switching-off for ? — > oo. Such a procedure has been 



introduced in [ ^3 1 as to describe intial correlations at some t = to evolving from an uncorrelated 
initial state at f — > — oo. In fact, the excess of the mass-shift contribution over its asymptotic value at 
finite times arises from spurious transient contributions that again lead to an unphysical UV scaling 
behavior and must be eliminated by switching the electromagnetic interaction off again according 



to (3.6) 



In addition to the vacuum contribution and the pure mass-shift contribution, the photon self- 



energy (3.10) also contains a contribution which describes the interference between the vacuum 



and the mass-shift effects. This contribution is also eliminated by the sequence of limits leading to 



(3.18). A more detailed consideration can be found in [38] 



6. Summary, conclusions and outlook 

In this work, we have investigated photon emission during the chiral phase transition in the 
early stage of a heavy-ion collision. During this phase transition, the quark mass is changed from 
its constituent value, ntc, to its bare value, mt, which leads to the spontaneous and non-perturbative 



pair production of quarks and antiquarks [ |39| , fK)| ]. This effectively contributes to the creation of 
the QGP, and we have investigated the photon emission arising from this pair-creation process. In 
particular, our investigations are relevant in the context of the question how the finite hfetime of 
the QGP during a heavy-ion collision affects the photon emission from it. 

Earlier investigations on this topic had been accompanied by two artifacts. In the first place, 
the photon number-density contained a divergent contribution from the vacuum polarization. Fur- 
thermore, the remaining contributions to this quantity were not integrable in the UV domain. In 
particular, the total number density and the total energy density of the emitted photons were loga- 
rithmically and linearly divergent, respectively. It has been our original motivation to resolve these 
conceptional problems in a satisfactory manner. 



For this purpose, we have first pursued a model approach [^ where the finite lifetime of 
the QGP is simulated by time-dependent quark/antiquai^k occupation numbers in the photon self- 
energy. By means of this procedure, we have been able to renormalize the divergent contribution 
from the vacuum polarization in a consistent manner, but it has not been possible to solve the 
remaining problem with the UV behavior to full extent. We have assumed that this remaining 
shortcoming results from a violation of the Ward-Takahashi identities. This aspect has been our 
eventual motivation to consider photon production arising from the change of the quark/antiquark 
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mass during the chiral phase transition since this scenario allows for a first-principle description by 
which the Ward-Takahashi identities ai^e conserved. 

In particular, the change of the quark/antiquark mass has been modeled by a Yukawa-like 
source term in the QED Lagrangian which couples the quarks and antiquarks to a scalar- and time- 



dependent background field. As in [BG, Bll, B3, 37], we have restricted ourselves to first-order QED 



processes. These processes are kinematically possible since the background field acts as a source 
of additional energy. On the other hand, the coupling to the background field has been resummed to 
all orders as to properly take into account the non-perturbative nature of the pair-creation process. 
In contrast to [3^ 31, 32, 37], we have not considered the photon number-density at finite times 



but for free asymptotic states. Such states have been obtained by switching the electromagnetic 
interaction according to the Gell-Mann and Low theorem. The photon number-density has then 
been considered in the asymptotic time limit ? — > oo and we have taken the switching parameter to 
zero at the very end of our calculation. By means of this procedure, the photon number-density 
does not contain any unphysical contributions from the vacuum polarization and is furthermore 
rendered UV integrable for physically reasonable mass parameterizations. 

In particular, we have shown that the consideration of this quantity for free asymptotic states 
is indeed cmcial to obtain such physically reasonable results. The reason is that a consistent def- 
inition actually is only possible for free asymptotic states, whereas a similar interpretation of the 
respective formal expression is usually not justified at finite times, t, since one then does not have 
free asymptotic states. The same problem occurs if the electromagnetic interaction is only switched 
on adiabatically from t -^ — oo but not switched off again for ? ^^ oo. This procedure has been sug- 
gested in p3| ] as to implement initial correlations at some t = to evolving from an uncorrelated 
initial state att ^- — oo. 

Moreover, the consideration of the 'photon number-density' at finite times effectively comes 
along with a violation of the Ward-Takahashi identities. Our investigations hence support the re- 



spective concern raised in [33, B4t] towards [30, 31 32] 



Consequently, our results indicate that the remaining artifact with the UV behavior within 



our model description [ |37| | results from an inconsistent definition of the photon number-density at 
finite times. This in turn imposes the question whether this artifact is also removed if the photon 
number-density is instead considered for free asymptotic states, which are obtained by switching 
the electromagnetic interaction according to the Gell-Mann and Low theorem, i.e., 

HEM{t)=efd^xY{x)Y^iV{x)A^'{x)^e-''\'\HEM{t), with £>0. (6.1) 



An accordant revision of []37[], for which the resulting photon spectra are depicted in Figure [11, 
shows that this is indeed the case. 



As in [ pQ , pl| , p2| ], the photon number-density decays as l/(Oj in the UV domain if the 



quai^k/antiquai^k occupation numbers are switched on instantaneously. If we turn from an instanta- 
neous switching to a switching over a finite interval of time, however, the photon number-density is 
suppressed from l/w^ to l/Wj which means that the total number density and the total energy den- 
sity of the radiated photons are both rendered UV finite. To the contrary, the UV scaling behavior 
is solely suppressed from ^/coS to l/coS'^ when turning from fi {t) to /2(0 if the 'photon number- 
density' is instead considered at finite times (see Figure ^. Our outlook to future investigations is 
hence as follows. 
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Figure 11: The photon spectra are also rendered UV integrable within our model description [^ if they are 
considered for free asymptotic states and if the time evolution of the QGP is modeled in a suitable manner, 
i.e., if its formation is assumed to take place over a finite interval of time. 



In the first place, the actual role of the Ward-Takahashi identities, which are violated in [ |370 
but consei"ved within our first-principle approach to chiral photon production requires deeper con- 
sideration. In particular, the fact that our model description ^l\\ leads to UV integrable photon 
spectra despite the violation of these identities does not necessarily disprove our earlier conjec- 
ture concerning their role. The reason is that the Ward-Takahashi identities can be violated in two 
different ways: 

• On the one hand, they can be violated directly by making ad-hoc assumptions on the two- 



time dependence of the photon self-energy, which has been the case in [37|. 



On the other hand, they can also be violated indirectly by considering the 'photon number- 
density' at finite times even though they are formally fulfilled at first. This has been shown 



in [38] 



Since our model description [ 37 1 leads to a UV integrable photon number-density if this quantity 
is considered for free asymptotic states, we hence have to determine if and, as the case may be, 
why only an indirect violation of the Ward-Takahashi identities, which occurs in addition to the 



direct violation when considering a transient 'photon number-density' within [|37|], leads to artificial 
results. 

Moreover, it is of particular interest whether our previous asymptotic description allows for a 
consistent extension to finite times and/or which alternative quantities can be considered to describe 
the time evolution of the electromagnetic sector during a heavy-ion collision. 
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